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Abstract. Nuclear inelastic scattering of synchrotron radiation has been used 
now since ten years as a tool for vibrational spectroscopy. This method has turned 
out especially useful in case of large molecules that contain a Mossbauer active 
metal center. Recent applications to iron-sulfur proteins, to iron(II) spin crossover 
complexes and to tin-DNA complexes are discussed. Special emphasis is given to 
the combination of nuclear inelastic scattering and density functional calculations. 
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1. Introduction 

Nuclear inelastic scattering (NIS) of synchrotron radiation (also known 
as nuclear resonant vibrational spectroscopy, NRVS) can be regarded 
as an extension of the conventional, energy-resolved Mossbauer spec- 
troscopy to energies of the order of molecular vibrations [1, 2]. NIS has 
similarity with resonance Raman spectroscopy: the energy of monochro- 
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matized synchrotron radiation is detuned from nuclear resonance by a 
certain amount of energy (several meV) that is sufficient to create or 
annihilate a phonon, or in terms of molecular spectroscopy, to excite 
or deexcite a molecular vibration. While the Raman intensity for a 
given molecular vibration depends on the change of polarizability with 
respect to the motion of the nuclei, the NIS intensity depends on 
the contribution of this particular vibration to the mean-square dis- 
placement (msd) of the Mossbauer nucleus. NIS is, therefore, a highly 
site-selective tool and it is as vibrational spectroscopy complementary 
to Raman or IR spectroscopy. The interpretation of measured NIS 
spectra is greatly facilitated due to the fact that the NIS intensity 
depends only on mechanical and not on electric properties of the inves- 
tigated molecules. For example, iron-ligand bond stretching vibrations 
of an iron-containing molecule are normally easy to identify because 
of their relatively large contribution to the msd. The dependence on 
purely mechanical properties also facilitates the simulation of measured 
NIS spectra with the help of quantum chemical calculations, which 
are usually based on density functional theory (DFT). The calculation 
of IR or Raman intensities, which depend on electric properties like 
changes of dipole moment or polarizability tensor, is far less accurate. 
A quantitative analysis of this fact with the example of the ferrocene 
complex is given in section 2. 

The comparison of experimental and simulated spectra provides the 
possibility to assign peaks in the measured NIS spectrum with calcu- 
lated normal modes of the molecule under study. This procedure has 
been applied first to small and medium-sized metal complexes (up to 
100 atoms) [3, 4] but is now applied also to larger biological molecules 
like rubredoxin [5] or DNA with chelated metal ions [6] (in the latter 
two cases DFT calculations are restricted to a cluster of about 100 
atoms around the iron center) . First applications of NIS to problems of 
biological relevance were related to model hemes like tetraphenylpor- 
phyrins [7, 8] and to heme-based proteins like myoglobin [7-9] (see also 
references in [10]). Recently NIS measurements have been performed 
on model complexes of iron-sulfur proteins [11], on a rubredoxin-type 
iron-sulfur protein in the reduced and in the oxidized state [5], and 
on Benzoyl-CoA reductase [5], which contains [4Fe-4S] 2+ centers. The 
NIS spectra for the rubredoxin-type protein allowed for the first time 
to investigate the change of the metal dynamics between the reduced 
and the oxidized form (see also sections 3.1 and 3.3). 

Another class of metal-containing molecules, which is being studied 
by NIS, is formed by the so called spin crossover (SCO) complexes. 
For these molecules NIS has turned out to be a very valuable method, 
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since the dynamics of the metal center governs the crossover from the 
low-spin to the high-spin isomer (section 3.4). 



2. Density functional theory calculations 

A detailed analysis of the dynamics of metal centers in large molecules 
requires combination of experimental investigations and theoretical cal- 
culations. In case of molecules that contain transition metals, den- 
sity functional theory (DFT) calculations are currently the method 
of choice. They have been shown to reproduce the frequencies and 
vibrational modes of transition metal compounds with reasonable accu- 
racy and can be applied also to large molecules containing hundreds of 
atoms. Frequencies and vibrational modes are used to simulate a NIS 
spectrum for comparison with experiment [3, 7]. Additionally, DFT 
calculations provide the partial vibrational density of states (PVDOS) 
for the Mossbauer nucleus which can be compared with a PVDOS that 
has been extracted from the measured spectrum. A detailed description 
how the calculated vibrational modes can be used to simulate NIS 
spectra can be found for instance in Refs. [3, 7]. 

In a low-temperature approximation the absorption probability S 
that is measured in a NIS experiment is proportional to the spectral 
contribution to the mean-square displacement (msd) of the Mossbauer 
nucleus: 

S(E,k) «^5(£-/^)((k- Ui ) 2 ) (1) 

i 

Here E and k denote energy and wavevector, respectively, of the ab- 
sorbed photon, index i enumerates calculated normal modes of vibra- 
tions, vi is the frequency of mode i and the vector Uj contains those 
three components of the normal coordinate vector that refer to the 
motion of the Mossbauer nucleus in x, y and z direction. The sum over 

all terms ^(k • Uj) 2 ^) yields the total measured msd of the Mossbauer 

nucleus. In order to simulate the NIS spectrum of a molecule it is, 
therefore, sufficient to know the frequencies V{ and to which degree 
the Mossbauer nucleus is participating in the individual vibrations. 
The latter information is given by the vectors Uj. The frequencies and 
normal modes of molecular vibrations may be regarded as mechanical 
properties of the molecule, although they are based on the molecule's 
electronic structure. 

With the example of the ferrocene molecule we demonstrate here 
that these mechanical properties are more easily calculated with suffi- 
cient accuracy than those electric properties (changes of dipole moment 
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Figure 1. NTS, IR and Raman intensities (solid, dotted and dashed line, respectively) 
of ferrocene in arbitrary units calculated with B3LYP/TZVP. 

or polarizability tensor) that are needed for the simulation of IR or 
Raman spectra. The dynamics of the iron center in ferrocene has been 
studied intensively by various spectroscopic techniques including NIS 
[12] and quasielastic nuclear forward scattering [13]. We have performed 
DFT calculations using the hybrid functional B3LYP [14] together 
with five different basis sets, 3-21G [15-20], CEP-31G [21-23], 6-311G 
[24, 25] (with the Wachters-Hay double zeta basis for Fe [26, 27]), TZV 
and TZVP [28, 29] as included in the Gaussian 03 program package [30]. 
The simulated spectra (Fig. 1) are in reasonable agreement with the 
experimental spectra given in ref. [12]. The comparison of calculated 
IR, Raman and NIS spectra in the region of metal-ligand vibrations 
exemplifies the complementarity of NIS and Raman spectra. In the 
region from to 60 meV nine vibrational modes are obtained from 
the calculation. Modes 1 and 7 are practically invisible for all three 
spectroscopic techniques, modes 2 and 3, modes 5 and 6 and modes 
8 and 9 are in pairs nearly degenerate. In this region mode 4 has the 
strongest Raman intensity, while the almost degenerate pair of modes 
8 and 9 exhibits the largest IR and NIS intensities. We have examined 
the influence of the basis set on the calculated IR, Raman and NIS 
intensities. In Fig. 2 the basis sets are ordered on the horizontal axis 
with increasing quality. It turns out that the calculated NIS intensity 
is almost independent of the employed basis set whereas the calcu- 
lated IR and Raman intensities exhibit large oscillations for small and 
medium sized basis sets. This result is even more striking if one takes 
into account that, according to Lambert-Beer law, there is an expo- 
nential relation between the calculated IR and Raman activities and 
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Figure 2. Relative NIS (•), IR (□) and Raman (A) intensities of mode 4 (Raman) 
and the nearly degenerate modes 8 and 9 (IR, NIS) of ferrocene (see text) calculated 
with B3LYP and different basis sets as indicated. For each spectroscopic technique 
the calculated intensities are divided by the value obtained with the TZVP basis 
set. The lines serve as guides to the eyes. 



the measured spectral intensities. In summary, a reasonable theoretical 
simulation of IR and, especially, of Raman spectra is computational 
much more demanding than a simulation of NIS spectra. This is an 
important aspect if large molecules are investigated where demanding 
computational methods are prohibitive. 



3. Recent examples 

3.1. RUBREDOXIN-TYPE IRON-SULFUR PROTEINS 

Iron-sulfur proteins are widely used for electron transport in cells. For 
a deeper understanding of the reorganisation of these proteins upon 
electron uptake or release the knowledge of the vibrational modes of the 
iron-sulfur core is of special importance. For the oxidized and reduced 
forms of a rubredoxin-type iron sulfur protein from Pyrococcus abyssi, 
which contains a single Fe-S4 core [31], we have measured NIS spectra 
at the nuclear resonance beamline ID 18 of ESRF in Grenoble, France. 
In the NIS spectrum of the oxidized form Fe-S bond-stretching modes 
could be identified which are in agreement with resonance Raman 
studies. The spectrum of the reduced form, where no Raman data are 
available, shows a down shift of these bond-stretching modes by about 
5 meV. This observation is in line with an increase of the Fe-S bond 
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length upon reduction as observed by EXAFS. In addition to the peaks 
that are assigned to the Fe-S bond stretching modes, a second group 
of prominent peaks was observed in the NIS spectra, which hardly 
change their positon upon reduction. In order to assign these peaks 
to molecular vibrations, DFT calculations have been performed for a 
series of models of the protein under study, Fe(SH)4, Fe(SCH3)4 and 
FeCys4, i.e. containing 9, 21 and 49 atoms, respectively. The simulated 
NIS spectra for these three models gain approach to the experimental 
spectra (oxidized and reduced) with growing model size and allow to 
assign the observed second group of peaks to Fe-S-C bending modes. 
The most simple model, Fe(SH)4, although far less realistic than the 
larger models, is nevertheless useful since it provides a qualitatively 
correct picture of the experimentally observed vibrations of the FeS4 
core. More details of this study are described in ref. [5]. 

3.2. Iron corroles 

Corroles are aromatic macrocycles that are closely related to porphyrins 
and have unique chemical properties such as stabilizing high oxidation 
states of coordinated metal ions, e.g. Fc(IV) [32]. We have performed 
NIS measurements on iron corroles with two different axial ligands at 
the nuclear resonance beamline ID 18 of ESRF. Currently, a detailed 
analysis of the observed molecular vibrations of these molecules is 
performed that exploits the complementarity of IR, Raman and NIS 
measurements in combination with DFT calculations [33] . 

3.3. TlN-DNA COMPLEXES 

NIS is a highly site-selective tool, focussed on Mossbauer nuclides. 
Therefore, most applications of NIS in biology and bioinorganic chem- 
istry have so far been performed with molecules containing iron centers. 
Recently, we have investigated tin-DNA complexes at the nuclear reso- 
nance beamline ID 22 of ESRF; this is, to our knowledge, the first ex- 
ample of a NIS study on a biological sample that contains a Mossbauer 
nuclide other than 57 Fe. This study has been performed to elucidate the 
binding site of organotin cations, Sn(II) and Sn(IV), chelated to DNA. 
The binding site of such tin-DNA complexes could not be determined by 
conventional structural methods like X-ray diffraction. Spectroscopic 
techniques like IR or Raman are hampered by the complexity of the 
spectra which are dominated by the vibrations of the large organic 
ligands. The site-selectivity of NIS focusses on those vibrations that 
are connected with a considerable msd of the tin nucleus. Our NIS 
investigation in combination with DFT calculations have corroborated 
the former hypothesis about the binding sites of tin that was based 
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on conventional Mossbauer studies and semiempirical molecular-orbital 
calculations [34]. More details can be found in refs. [6, 35]. 

3.4. Spin crossover complexes 

These molecules exhibit a transition from a low-spin (LS) to a high- 
spin (HS) isomer that is induced by change of temperature or pressure 
or by irradiation with light [36] . Spin crossover complexes are not only 
very promising materials for technical applications, they also serve as 
models for the transition of the metal spin in biomolecules, observed for 
instance in a catalytic cycle. The driving force for the spin transition 
is the entropy difference between the two spin isomers. The major 
contribution to this entropy difference arises from the dynamics of 
the metal ion: in the LS state the d electrons of the metal center 
occupy non-bonding or weakly 7r-anti-bonding molecular orbitals. Upon 
transition to the HS isomer part of the d electrons are transferred to 
(T-anti-bonding molecular orbitals. Consequently, in the HS isomer the 
metal-ligand bonds are weaker and longer and also softer, i.e. the force 
constants of the metal-ligand bond stretching vibrations are lower. This 
leads to a significant down-shift of the corresponding frequencies upon 
crossover from LS to the HS state (about 10 to 20 meV in case of 
typical Fe(II) SCO complexes). For this reason, the HS isomer has - in 
comparison to the LS isomer - many more vibrational states that are 
thermally accessible, and the vibrational entropy contribution of the HS 
isomer is significantly higher than the corresponding value of the LS 
isomer (some 10 J mol -1 K _1 in case of Fe(II) SCO complexes). First 
NIS experiments have been performed on samples of [Fe(tpa)(NCS)2] 
(tpa = tris(2-pyridylmethyl)amine) and [Fe(bpp)2](BF4)2 (bpp = 2,6- 
bis(pyrazol-3-yl)pyridine) [3, 37, 38]. Subsequent measurements on a 
single crystalline sample of [Fe(tptMetame)] (0104)2 (tptMetame = 1,1,1- 
tris((N-(2-pyridylmethyl)-N-methylamino)methyl)ethane) yielded an- 
gular resolved NIS spectra [4]. With the same sample the influence 
of a pressure induced SCO on the metal dynamics has been investigated 
with NIS [39]. Studying [Fe(phen) 2 (NCS) 2 ] (phen = 1,10-Phenanthroline) 
the complementarity of NIS and IR and Raman spectroscopy has been 
demonstrated. Combining these three vibrational spectroscopic tech- 
niques with DFT calculations allows an almost complete assignment 
of the observed vibrational modes [40]. One of the most recent NIS 
experiments on SCO complexes [41] has been used to investigate five 
different phases of the complex [Fe(ptz)e](BF4)2 (ptz = 1-n-Propyl- 
tetrazole) that can be reached by combination of temperature changes 
and irradiation with light of different frequencies [42] . 
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4. Conclusion 

In summary we conclude that the comparison of complementary ex- 
perimental (i.e. NIS, IR, Raman) and theoretical spectra provides the 
possibility to assign measured peaks with calculated normal modes of 
the metal complex under study. This procedure has been applied to 
small and medium-sized iron complexes and has recently been applied 
also to iron-containing proteins and to DNA being chelated with tin. 
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